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Abstract 

A  semi-empirical  equation  was  proposed  to  simulate  the  behavior  of  a  direct  methanol  fuel  cell  (DMFC).  The  individual  voltage  losses  in 
a  DMFC  due  to  methanol  crossover  and  the  overpotentials  of  both  the  cathode  and  anode  can  be  distinguished.  Three  sets  of  experiments  were 
designed  and  carried  out  to  account  for  the  three  voltage  losses.  The  values  of  each  parameter  in  the  model  were  then  calculated  and  the  computation 
showed  that  the  fitted  result  and  the  experimental  data  were  well  matched.  The  relation  between  each  significant  phenomena  and  each  parameter 
is  discussed.  The  model  quantitatively  identified  the  major  voltage  losses  to  be  both  the  sluggish  reaction  of  methanol  oxidation  on  the  anode  and 
the  slow  oxygen  reduction  on  the  cathode.  The  impact  on  cell  performance  by  manipulating  individual  parameters  is  also  discussed. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  DMFC  is  a  good  potential  portable  power  source  due  to 
its  low  operating  temperature  (room  —  60  °C)  and  high  theoret¬ 
ical  energy  density  (~5000  WhL-1).  At  present,  a  low  power 
density  and  a  high  catalyst  loading  of  the  membrane  electrode 
assembly  (MEA)  are  the  major  barriers  that  inhibit  its  commer¬ 
cialization.  To  improve  the  MEA  power  density,  we  need  to 
analyze  the  MEA  behavior  and  identify  the  critical  factors  that 
limit  power  density.  However,  the  behavior  inside  the  MEA  is 
extremely  complicated  and  heavily  dependent  on  the  operating 
condition,  the  constituting  materials,  the  MEA  structure,  and  the 
manufacturing  procedure,  etc.  [1-6].  Thus,  theoretical  models 
or  semi-empirical  equations  have  been  developed  in  order  to 
predict  and  analyze  factors  which  affect  the  performance. 

The  power  curve  of  proton  exchange  membrane  fuel  cells 
(PEMFC)  have  been  extensively  studied.  In  general,  they  can 
be  divided  into  three  regions,  the  active,  the  ohmic  and  the 
mass-transfer  limiting  regions  [7-10].  In  the  low  current  den¬ 
sity  region,  the  activation  overpotential  of  oxygen  reduction  at 


*  Corresponding  author.  Tel.:  +886  3  5715131  3660;  fax:  +886  3  5715408. 
E-mail  address:  ccwan@mx.nthu.edu.tw  (C.-C.  Wan). 

0378-7753/$  -  see  front  matter  ©  2006  Elsevier  B.V.  All  rights  reserved, 
doi :  10. 1016/j .  j  po  wsour. 2005. 12.052 


the  cathode  predominates.  In  the  intermediate  region,  the  cell 
internal  resistance,  mainly  contributed  by  the  membrane  resis¬ 
tance,  becomes  the  major  factor  resulting  in  a  linearly  decreasing 
section  of  the  discharge  curve.  In  the  high  current  density  region, 
the  overall  cell  reaction  rate  is  limited  by  the  depletion  of  reac¬ 
tants.  Consequently,  the  mass-transfer  overpotential  becomes 
the  dominant  factor  and  it  induces  a  sharp  decay  in  the  power 
density. 

Many  semi-empirical  models  have  been  proposed  to  describe 
the  fuel  cell  power  curve  [2,3,7-10].  These  models  serve  to 
distinguish  the  potential  drop  terms  individually.  Table  1  sum¬ 
marizes  several  important  semi-empirical  equations  for  the 
PEMFC.  In  1988,  Srinivasan  et  al.  [7]  presented  a  semi-empirical 
equation  to  describe  the  discharge  behavior  of  a  PEMFC.  This 
equation  contains  a  logarithmic  term  for  the  activation  over¬ 
potential  and  a  linear  term  for  the  ohmic  drop.  This  equation 
is  based  on  the  Tafel  equation  and  can  perfectly  describe  the 
polarization  curve  before  reaching  the  mass-transfer  limiting 
region.  Kim  et  al.  [8]  then  added  an  exponential  term  to  account 
for  the  mass-transfer  overpotential.  Although  Kim’s  equation 
can  fit  the  whole  curve  well,  this  exponential  term  is  merely 
data  fitting  and  has  no  physical  meaning.  So  it  was  revised  by 
Squadrito  et  al  [9]  and  Chu  and  Jiang  [10]  with  some  theoretical 
reasoning. 
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Nomenclature 

Symbols 

^3-cat 

Actual  available  catalyst  area 

^MEA 

The  MEA  cross-section  area 

B 

Tafel  slope 

Gin,ch 

The  methanol  concentration  on  the  anodic  GDL 
surface 

Gin, cat 

The  methanol  concentration  at  the  interface 
between  GDL  and  the  catalyst  layer 

cef 

The  reference  methanol  concentration 

E 

Electrode  potential 

Gcell 

Cell  voltage 

Across 

The  intrinsic  voltage  loss  due  to  methanol 

crossover 

£r,ME 

The  reversible  electrode  potential  of  methanol 
oxidation 

£r,02 

The  reversible  potential  of  oxygen 

Eo 

Defined  in  Eq.  (7)  for  anodic  part  and  Eq.  (13)  for 
cathodic  part 

E 

Faradaic  constant 

I 

Current 

i 

The  apparent  current  density 

i* 

Lo 

The  apparent  exchange  current  density  at  the  ref¬ 
erence  concentration  of  methanol  or  partial  pres¬ 
sure  of  the  oxygen 

*  cross  ocp  The  aPParent  current  density  resulting  from 
crossover  at  the  open  circuit  potential 

&cl,ME 

The  methanol  mass-transfer  coefficients  inside 
GDL 

kcff 

The  effective  mass-transfer  coefficient  inside  the 
GDL 

m 

Equal  to  the  reaction  order  multiplied  by  the  Tafel 
slope 

N 

The  mass  flux  of  methanol 

n 

The  inverse  of  the  apparent  limiting  current  den¬ 
sity 

R  in,ch 

The  oxygen  partial  pressure  at  the  GDL  surface 

R  in, cat 

The  partial  pressure  at  the  interface  between  GDL 
and  cathodic  catalyst  layer 

pref 

The  reference  partial  pressure  of  oxygen 

^cell 

The  internal  cell  resistance 

r 

The  reaction  order 

z 

The  equivalent  of  specie 

Greek  symbols 

a 

The  charge-transfer  coefficient 

8 

The  proportional  constant 

0 

An  empirical  factor  which  represents  the  ratio  of 
the  oxygen  reversible  potentials  with  and  without 
methanol  existing 

Suffix 

a 

Anode 

c 

Cathode 

Table  1 

The  empirical  equations  for  a  PEMFC 

Equation  References 

E  =  Eo  —  b  log  i  —  Ri  Srinivasan  et  al.  [7] 

E  =  Eo  —  b  log  i  —  Ri  —  m  exp (ni)  Kim  et  al.  [8] 

E  =  Eo  —  b  log  i  —  Ri  —  odk  ln(l  —  f>i)  Squadrito  et  al.  [9] 

E  =  Eq  —  b  log  i  —  Ri  —  imm  exp (nim)  Chu  and  Jiang  [10] 


The  same  approach  has  been  applied  to  the  DMFC.  Scott 
et  al.  [11]  proposed  a  simplified  equation  in  which  the  anodic 
mass-transfer  overpotential  was  represented  by  the  concentra¬ 
tion  overpotential  of  methanol.  Scott  and  coworkers  [12]  found 
that  Kim’s  equation  was  unsatisfactory  in  the  low  concentration 
case.  In  contrast,  Squadrito ’s  equation  did  show  perfect  fittings 
under  all  testing  conditions  but  it  needed  heavy  computing  to 
solve  the  equations.  Thus,  Scott’s  group  presented  a  modified 
equation  as  follows: 

£ceii  =  £o  -  b  log  i  -  Rei  +  Ci  ln(l  -  C2i)  (1) 

In  2003,  Scott  and  coworkers  [13]  based  on  the  Tafel  approx¬ 
imation  developed  a  formula  similar  to  Eq.  (1),  in  which  each 
parameter  (Eq,  b ,  Re,  C\  and  C2)  has  its  physical  meaning.  Eval¬ 
uation  of  individual  parameters  could  help  identify  the  key  factor 
that  causes  decay  of  cell  performance. 

However,  most  semi-empirical  equations  could  not  separate 
the  anodic  voltage  loss  from  the  cathodic  voltage  loss.  This 
is  acceptable  for  the  PEMFC  because  the  anodic  overpoten¬ 
tial  is  minor  and  can  be  ignored.  Unlike  the  PEMFC,  methanol 
crossover  from  anode  to  cathode  and  the  sluggish  reaction  of 
methanol  oxidation  on  the  anode  are  two  unique  characteristics 
of  a  DMFC.  The  methanol  crossover  is  due  to  the  diffusion  and 
the  osmostic  drag  of  methanol  from  anode  to  cathode.  It  not 
only  lowers  the  output  voltage  of  a  DMFC  but  also  reduces  the 
utilization  efficiency  of  methanol  fuel.  The  sluggish  oxidation 
of  methanol  causes  an  additional  activation  overpotential  loss  at 
the  anode.  Consequently,  the  equations  listed  in  Table  1  are  not 
adequate  for  a  DMFC. 

Thus,  we  have  developed  a  semi-empirical  model  for  a  DMFC 
by  taking  these  two  important  factors  into  account.  Moreover, 
this  model  could  separate  the  anodic  and  the  cathodic  perfor¬ 
mances  individually  with  a  clear  crossover  effect.  Major  voltage 
losses  of  the  DMFC  can  then  be  identified  from  its  power  deliv¬ 
ery  curve. 

2.  Theoretical  derivation 

Fig.  1(a)  is  a  basic  sketch  of  the  MEA  structure  of  a  DMFC 
including  flow  fields,  gas  diffusion  layers  (GDLs),  catalyst  lay¬ 
ers  on  both  electrodes  and  the  membrane.  The  GDL  is  made  of 
carbon  paper  or  carbon  cloths  that  have  large  pores  for  fast  dif¬ 
fusion  of  species.  The  catalyst  layers  have  a  fine  pore  structure 
and  a  large  surface  area  so  that  the  species  react  with  a  min¬ 
imum  amount  of  catalyst.  At  the  anode-membrane  interface, 
the  residue  methanol  will  cross  to  the  cathode  due  to  osmostic 
drag  and  diffusion.  Fig.  1(b)  depicts  the  methanol  concentration 
distribution  on  the  anodic  side  and  the  oxygen  pressure  on  the 
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Fig.  1.  (a)  A  simplified  structure  of  the  DMFC  MEA;  (b)  the  methanol  concen¬ 
tration  and  oxygen  pressure  profile  at  different  positions. 


cathodic  side.  At  steady  state,  two  assumptions  were  made  inside 
the  GDL:  (i)  only  the  mass  transport  of  species  occur  and  (ii) 
the  reactant  flux  is  proportional  to  the  concentration  gradient. 

Inside  the  anodic  GDL,  the  mass  flux  of  methanol  (Ame)  is 
directly  proportional  to  the  difference  of  methanol  concentration 
on  the  GDL  surface  (Cin,ch)  and  at  the  interface  between  GDL 
and  the  catalyst  layer  (Cin,cat). 

A^ME  —  ^cl,ME(Gin,ch  Gin,cat)  (2) 

where  &ci,me  is  the  methanol  mass-transfer  coefficients  inside 
GDL. 

If  methanol  crossover  flux  does  not  exist,  all  methanol  which 
diffuses  from  the  GDL  should  be  electrochemically  oxidized  to 
carbon  dioxide  on  the  anodic  catalyst  layer.  However,  a  DMFC 
has  significant  methanol  crossover  from  the  anode  to  the  cathode 
[2-4, 14-17] .  The  amount  of  methanol  crossover  could  be  treated 
as  linearly  decreasing  with  increasing  current  density  [18-21]. 
The  material  balance  of  methanol  on  the  anode  can  be  described 
by  Eq.  (3).  The  amount  of  methanol  diffusing  from  the  GDL 
to  the  anodic  catalyst  layer  (AmeAmea)  is  equal  to  the  sum  of 

(i)  methanol  consumed  in  the  catalyst  layer  (7a  in  Eq.  (3))  and 

(ii)  the  amount  of  methanol  diffusing  and  migrating  to  cathode 

(Across, OCP  —  ^a). 

XT  A  4f  (Across, OCP  —  <^a)  /ox 

NmeAmea  = - — -  (3) 

z&F 

where  /cross, ocp  is  the  equivalent  current  due  to  crossover  at  the 
open  circuit  potential,  8  a  ratio  describing  how  the  crossover  flux 
decreases  with  increasing  current,  za  the  chemical  equivalent  of 
methanol,  F  Faradic  constant,  and  Amea  the  cross-section  area 
of  MEA,  respectively. 

The  anodic  current  density  should  be  equal  to  the  anodic  cur¬ 
rent  (/a)  divided  by  the  actual  catalyst  area  (Acat,a).  However,  the 
total  active  area  is  difficult  to  measure  and  it  changes  with  various 
manufacturing  methods.  Therefore,  the  MEA  cross-section  area 
(Amea)  was  used  in  our  analysis  instead  of  Acat,a.  Accordingly, 
the  current  densities  in  Eq.  (4)  become  the  apparent  current  den¬ 


sities  (/*  and  /*ross  OCP),  which  are  equal  to  the  current  divided 
by  Amea- 


Nme  = 


it  +  d 


cross,  OCP 


-  K) 


zaF 


Substituting  Ame  from  Eq.  (2)  into  (4),  Eq.  (5)  was  obtained. 


(1  5)/a  +  *  cross,  OCP  —  ^cl,ME£aG(Cin,ch  Gin,cat) 


The  apparent  current  density  (/*)  and  anodic  potential  (Ga)  is 
assumed  to  follow  the  Tafel  behavior  as  given  by  Eq.  (6): 


•  >k  .;k 

7  =  l 

loa 


c 


in, cat 


cref 


exp 


^zaF 

RT 


(Ga  —  Gr?  me) 


where  Cref  is  the  methanol  concentration  at  a  fixed  reference 
condition,  Gr, me  the  reversible  electrode  potential  of  methanol 
oxidation,  ra  the  reaction  order  of  methanol,  aa  the  charge- 
transfer  coefficient  of  methanol  oxidation,  and  /*a,  the  exchange 
current  density  of  methanol  oxidation  at  a  reference  methanol 
concentration,  C 7ef . 

Eq.  (5)  coupled  with  Eq.  (6)  can  be  rearranged  into  Eq.  (7), 
in  which  Ga  is  expressed  as  a  function  of  /*. 

Ga  =  G0, a  +  ba In  i*  -  maln{l  -  wa[(  1  -  <$)/*  +  Coss,ocp]} 

(7) 


where 

1 

Go, a  =  Gr,ME  +  In  t— 

7 * 

loa 


rref 


RT 

aazaF’ 


raRT  1  1 

- >  na  ~  -  =  - 

WaZaF  ^aGfeci,MEGin,ch  him,  a 


The  variable  8  is  defined  to  be  the  ratio  of  /*ross  OCP  to  q*m  a  (or 
1  /na)  as  given  by  Eq.  (8).  The  methanol  crossover  flux  has  the 
largest  amount  at  the  open  circuit  condition  (/*ross  ocp)’  ^ut  ^ 
approaches  0  when  /*  reaches  the  methanol  mass-transfer  lim¬ 
iting  current  density  (q*m  a)  and  Cin,cat  approaches  0. 

^  ^a^cross,OCP  (8) 


When  we  combine  Eq.  (8)  with  Eq.  (7),  the  anodic  potential  of 
DMFC  can  be  expressed  as  follows: 


Ea  =  £o,a  -  Across, a  +  ba  In  i*  -  ma  ln(l  -  na/*)  (9) 


where 


Gcross,a  —  Wta  ln(l  8) 

The  cathodic  current  density  (/*)  can  also  be  obtained  sim¬ 
ilarly.  The  material  balance  of  oxygen  on  the  cathode  can  be 
described  by  Eq.  (10)  where  £eff,c  is  the  effective  mass-transfer 
coefficient  of  oxygen  inside  the  cathodic  GDL  and  zc  is  the  chem¬ 
ical  equivalent  of  oxygen.  The  right-hand-side  of  Eq.  (10)  repre¬ 
sents  the  oxygen  flux  in  GDL  and  is  proportional  to  the  difference 
of  oxygen  partial  pressure  at  the  GDL  surface  (Pin,ch)  and  the  par¬ 
tial  pressure  at  the  interface  between  GDL  and  cathodic  catalyst 
layer  (Gin,cat).  We  visualized  that  the  oxygen  transported  from 
GDL  is  reduced  to  /*  and  consumed  by  the  oxygen/methanol 
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redox  reaction  in  the  form  of  a  local  cell  which  should  equal  to 
/*ross  ocp  —  Si * .  So  a  material  balance  of  oxygen  leads  to  Eq. 


(10). 


ic  ”1“  Across, OCP  Sic)  —  keff,cZcF(Pm,ch  ^in,cat)  (10) 


The  oxygen  reduction  proceeds  in  parallel  on  the  methanol- 
covered  and  fresh  catalyst  surfaces  if  methanol  passes  through 
the  membrane  [22] .  On  the  methanol-free  surface,  oxygen  reacts 
with  the  protons  and  electrons  are  released  at  the  anode,  as  shown 
in  Eq.  (11). 


•  >k  •  >k 

1  =  1 
%  Loc 


rin,cat 

pvQf 


exp 


OicZcF 

RT 


( Eq  ErQ2) 


(11) 


where 


rcRT 

nc 

OicZcF 


1 


Zc  FJCq  ffjC  P\n,ch 


RT 

OicZcF' 


The  apparent  current  density  flowing  out  of  the  cathode  (/*)  is 
equal  to  the  apparent  current  density  flowing  out  of  the  anode 
(/*).  To  simplify  the  notation,  we  used  the  apparent  current  den¬ 
sity  f  for  i*  in  Eq.  (13)  and  for  /*  in  Eq.  (9).  The  cell  voltage 
(Eceii)  is  the  difference  between  Ec  and  Ea  if  there  is  no  inter¬ 
nal  resistance  (ftceii).  When  RCc\\  is  taken  into  account,  the  cell 
voltage  Eceii  is  as  follows: 


where  Eref  is  the  oxygen  partial  pressure  at  reference  condi¬ 
tion.  Ey  o2  is  the  reversible  potential  of  oxygen.  ac  and  rc  are 
the  charge-transfer  coefficient  and  the  reaction  order  of  oxygen 
reduction,  respectively,  and  /*c  is  the  apparent  exchange  current 
density  of  oxygen  reduction  at  Eref . 

On  the  methanol-covered  surface,  the  methanol  oxidation  and 
oxygen  reduction  form  a  redox  couple.  The  amount  of  methanol 
crossover  from  the  anode  (/*ross  OCP  —  8i*)  is  totally  consumed 
at  the  cathode.  Since  the  methanol  is  oxidized  at  much  lower 
potential  than  the  potential  of  oxygen  reduction,  whereby  form¬ 
ing  a  local  cell  reaction.  The  presence  of  methanol  accelerates  the 
oxygen  reduction  and  it  should  have  a  lower  specific  free  energy 
on  the  methanol-covered  catalyst  surface  [22] .  On  this  surface, 
the  methanol  oxidation  rate  is  equal  to  the  oxygen  reduction  rate 
as  given  by  Eq.  (12). 


•  r\  • 

*  cross,  OCP  °lc  loc 


P  \rc 
rm,cat  \ 

pref  J 


x  exp 


OicZcF 


RT 


0ETio2) 


(12) 


where  the  empirical  factor,  6 ,  is  always  smaller  than  1.  Obvi¬ 
ously,  when  6  equals  unity,  Eq.  (12)  returns  to  Eq.  (11).  It 
represents  the  diminution  of  free  energy  from  a  methanol-free 
surface  to  a  methanol-covered  surface. 

Substituting  i*  from  Eq.  (11)  and  z*ross  OCP  —  <$/*  from  Eq. 
(12)  into  Eq.  (10),  Eq.  (13)  is  obtained: 

Ec  =  £o,c  +  bc  In  1 1  -  exp 

—  be  ln[(l  -  8)i*  +  * *ross, OCP ] 

+  mc  ln{  1  -  nc[(  1  -  <$)/*  +  /*ross?0cp]}  (13) 


(1  ~  0)Er,O2 
be 


FcqW  —  Fc  Ea  ^celfl 

—  (Eo,c  T  Across, c)  —  be  In [( 1  —  8)i  +  ^cross,OCp] 

+  rnc  ln{  1  -  nc[(  1  -  5)/*  +  /*ross,OCp]}  “  E0,SL 

T  Across, a  1^(^  )  H-  ^a  ln(l  Jl^i  )  ^celfl 

(14) 

where  we  define  a  term  £cross?c,  which  equals  ln{  1  —  exp[— (1  — 
6)EY^o2/bc]},  and  its  physical  meaning  will  be  explained  later. 

3.  Experimental  and  numerical  fitting  procedure 

The  MEA  used  in  our  experiment  had  a  carbon  cloth  as  the 
GDL  on  both  sides  and  its  size  was  10.89  cm2  (3.3  cm  x  3.3  cm). 
The  anodic  catalyst  layer  had  4  mg  cm-2  PtRu  black  with  a  1 : 1 
atomic  ratio.  The  cathodic  catalyst  layer  contained  4  mg  cm-2 
Pt  black.  The  proton  exchange  membrane  was  Nation  117  from 
DuPont. 

In  order  to  calculate  each  parameter  in  Eq.  (14),  we  did 
three  sets  of  experiments.  Experiment  A  (Exp  A)  was  actually  a 
PEMFC  system  fed  with  hydrogen  at  the  anode.  Since  the  anodic 
polarization  can  be  neglected  [7-10],  Exp  A  can  in  fact  give  a 
measure  of  the  cathodic  performance  of  DMFC  under  the  same 
condition.  Experiment  B  (Exp  B)  provides  regular  DMFC  dis¬ 
charging  with  methanol  and  air  as  feed.  In  Experiment  C  (Exp 
C),  the  cathode  was  fed  with  hydrogen  so  that  cathodic  polar¬ 
ization  can  be  neglected  [23]  and  the  anodic  performance  of  the 
DMFC  can  be  identified.  Table  2  listed  the  conditions  used  in 
each  experiment.  We  chose  flow  rates  higher  than  conventional 
figures  in  order  to  keep  the  input  concentration  constant  and  our 
assumption  in  the  theoretical  derivation  that  the  reactant  transfer 
is  only  dependent  on  diffusion  and  independent  of  the  input  flow 
rate  is  thereby  justified. 


Table  2 


The  flow  rate  and  cell  temperature  settings  in  the  experiments  (10%,  v/v  ~  2.5  M) 


Experiment 

Anode 

Cathode 

Cell  temp.  (°C) 

Reactant 

Flow  rate  (1  min  1 ) 

Reactant 

Flow  rate  (1  min  1 ) 

Exp  A 

h2 

1 

Air 

1 

40,60 

Exp  B 

CH3OH  (10%,  v/v) 

0.15 

Air 

1 

40,60 

Exp  C 

CH3OH  (10%,  v/v) 

0.15 

h2 

1 

40,60 
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The  current  scan  rates  and  ranges  for  the  different  experiments 
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Experiment 

Method 

Low  current 

Large  current 

Scan  rate  (mA  s“ 

_1)  Scan  range  (A) 

Scan  rate 

Scan  range  (A) 

Exp  A 

0.1 

0-0.3 

0. 1  A  min- 1 

0.3-10 

Exp  B 

0.1 

0-0.3 

0.025  A/30  s 

0.3-3.5 

Exp  C 

0.1 

0-0.3 

1  mAs-1 

0.3-0.7 

AC  impedance  (at  OCP) 

To  measure  the  cell  resistance  individually  under  open  circuit  potential  after  every  test 

In  these  three  experiments,  the  measuring  instruments  and 
conditions  in  the  low  current  region  were  different  from  those 
in  the  high  current  density  region.  Table  3  lists  the  measuring 
conditions  used  in  these  two  current  ranges.  In  Exp  A  and 
Exp  B,  a  potentiostat  (Solartron  1280z)  was  employed  to 
obtain  the  precise  variation  in  the  activation  region  at  low 
currents  (0-0.3  A)  and  the  fuel  cell  test  station  (Electrochem. 


Fuelcell  2000)  provided  the  polarization  curve  in  the  high 
current  region  (above  0.3  A).  However,  Exp  C  was  totally 
done  by  the  potentiostat  because  of  the  limitation  of  the 
test  station.  We  also  measured  the  cell  resistance  of  each 
experiment  at  open  circuit  by  AC  impedance.  The  internal 
resistance  of  Exp  B  was  used  as  the  cell  resistance  (RCq\i)  in  our 
analysis. 


Fig.  2.  The  calculation  procedure  of  the  semi-empirical  equation.  R\,  Rb  and  Rq  are  the  cell  resistances  of  Exp  A,  Exp  B  and  Exp  C,  respectively.  Moreover,  we 
also  use  Rb  as  Rce n  in  our  semi-empirical  model. 
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The  parameters  in  Eq.  (14)  could  be  calculated  by  using 
least-square  curve  fitting.  Cathodic  parameters  (£b,c>  bc,  Ecr0ss, c, 

'cross.OCP-  mc  and  nc)  and  anodic  parameters  (£0,a.  ba,  Ec ross,a, 
raa  and  fta)  were  calculated  in  a  sequential  manner  from  these 
three  sets  of  data  mentioned  above.  The  calculation  procedure 
is  shown  in  Fig.  2.  The  internal  resistance  of  Exp  A  ( R\ ),  Exp  B 
( Rb ),  and  Exp  C  (Rq)  were  measured  by  the  AC  impedance 
technique.  The  ohmic  overpotential  resulting  from  the  inter¬ 
nal  resistance  was  deducted  before  least- square  curve  fitting. 
Because  the  ohmic  drop  is  beforehand  eliminated,  we  only 
need  to  consider  the  activation  overpotential  and  the  concen¬ 
tration  overpotential  in  our  calculation.  Accordingly,  the  effect 
of  the  concentration  overpotential  is  enhanced  in  the  polariza¬ 
tion  curve,  so  we  can  evaluate  m  and  n  even  from  data  not  totally 
in  the  limiting  region. 

Firstly,  we  calculated  the  cathodic  parameters,  £o,c>  &c> 
and  nc  from  data  of  Exp  A.  Due  to  the  low  current  range  of  Exp 
C,  (£o,a  —  Ecr0Ss,a)  and  Z?a  could  be  obtained  by  neglecting  the 
concentration  overpotential  or  ma  and  rza.  The  sum  of  Exp  B’s 
and  Exp  C’s  polarization  curves  should  be  the  actual  cathodic 
performance  of  DMFC.  In  contrast,  the  polarization  curve  of 
Exp  A  should  be  the  ideal  cathodic  performance  because  it  has 
no  crossover  effect.  Therefore,  the  cathodic  parameter  related 
to  methanol  crossover,  £cr0Ss,c  and  i *ross  QcP’  could  be  evalu¬ 
ated  by  fitting  the  deviation  between  ideal  and  actual  cathodic 
polarization  curve.  Finally,  the  difference  between  the  potential 
response  calculated  with  the  parameters  already  obtained  and 
actual  data  of  Exp  B  is  represented  by  the  maln(l  —  n^i  )  term  in 
Eq.  (14).  So  ma  and  n a  can  be  evaluated  from  this  discrepancy. 
Moreover,  £o,a  and  Ec r0Ss,a  are  also  separated  easily  according 
to  the  definition  of  Ecr0Ss,a  shown  in  Eq.  (9). 

4.  Results  and  discussions 

4.1.  Model  validations 

The  measured  data  of  cell  voltage  versus  current  density  (£ceii 
versus  /*)  for  different  schemes  are  given  in  Figs.  3-5.  The  val- 


Fig.  4.  The  curve  fitting  result  of  Exp  B. 


ues  of  the  parameters  listed  in  Tables  4  and  5  were  obtained 
by  fitting  Eq.  (14)  with  the  curves  in  the  above  figures  except 
the  resistances.  The  £ceii  versus  C  curves  calculated  from  these 
parameters  are  also  plotted  on  the  same  diagram  for  compari¬ 
son.  This  semi-empirical  model  obviously  fits  the  experimental 
data  very  closely  with  H2/air  feeding,  methanol/air  feeding,  or 
methanol/H2  feeding  at  40  and  60  °C.  It  should  be  noticed  that 

Table  4 

Parameters  of  semi-empirical  model  obtained  by  fitting  the  data  at  40  °C 


Eo, c  (V) 

0.817 

bc  (Ve-1) 

0.021 

mc  (V  e-1) 

0.074 

nc  (cm2  A- 1 ) 

1.020 

Ecross,c  (V) 

-0.065 

Cross, OCP  (A  Cm  2) 

0.030 

£b,a  (V) 

0.493 

ba  (Ve-1) 

0.031 

raa  (Ve-1) 

0.051 

na  (cm2  A- 1 ) 

2.660 

-Across, a  (V) 

-0.004 

R  (Q  cm2) 

0.563 

Fig.  3.  The  curve  fitting  result  of  Exp  A. 
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Table  5 

Parameters  of  semi-empirical  model  obtained  by  fitting  the  data  at  60  °C 


£o,c  (V) 

0.827 

bc  (Ve-1) 

0.021 

mc  (V  e-1) 

0.056 

nc  (cm2  A- 1 ) 

0.997 

-Across, c  (V) 

-0.064 

Cross, OCP  ^  cm  2) 

0.114 

£b,a  (V) 

0.406 

A  (Ve-1) 

0.030 

raa  (Ve-1) 

0.131 

na  (cm2  A- 1 ) 

1.607 

-Across, a  (V) 

-0.027 

R  (Q  cm2) 

0.579 

the  apparent  resistance  in  Fig.  3  (slope  in  the  linear  region)  was 
obtained  with  the  membrane  saturated  with  water  (for  hydro¬ 
gen/air  feeding).  The  resistances  listed  in  Tables  4  and  5  were 
measured  by  the  impedance  method  with  the  membrane  sat¬ 
urated  with  a  mixture  of  methanol  and  water  (methanol/air 
feeding)  due  to  methanol  crossover.  The  resistance  of  the  mem¬ 
brane  saturated  with  water  was  lower  than  the  resistance  of  the 
membrane  saturated  with  methanol/water.  This  is  reflected  in  the 
resistance  difference  between  Fig.  3  and  Tables  4  and  5.  Several 
observations  were  made  as  follows: 

(i)  Because  a  higher  temperature  causes  faster  reaction  kinet¬ 
ics  of  both  methanol  oxidation  and  oxygen  reduction,  the 
activation  overpotential  should  decrease  with  increasing 
temperature.  In  other  words,  a  lower  £o,a  and  a  higher 
Eq  c  are  obtained  when  the  temperature  rises.  Fitting  results 
demonstrate  the  same  trend  as  shown  in  Tables  4  and  5. 

(ii)  In  the  literature  [4,16,22],  Z?a  for  methanol  oxidation  is 
in  the  range  of  60-120  mV  dec-1  (25.6-51.2  mV  e-1)  and 
bc  for  oxygen  reduction  should  be  around  60  mV  dec-1 
(25.6  mV  e-1).  As  shown  in  Tables  4  and  5,  the  value  of 
Z?a  is  in  the  predicted  range,  but  the  value  of  bc  is  slightly 
lower.  It  is  noteworthy  that  Z?a  and  bc  are  insensitive  to  the 
temperature  variation  within  our  experimental  range.  (The 
difference  of  Z?a  or  bc  between  40  and  60  °C  is  only  about 

2  mV  e-1.)  Therefore,  the  polarization  is  sensitive  to  tem¬ 
perature  due  to  the  changes  of  £o,a  and  Across, a  (or  ^o,c  and 
Across, c  at  the  cathode)  at  different  thermal  settings.  This 
will  be  discussed  later. 

(iii)  ma  and  mc  are  related  to  mass  transfer  and  are  defined  to 
be  the  reaction  order  multiplied  by  the  Tafel  slope.  In  the 
literature  [18,19,24-26],  the  reaction  orders  of  oxygen  and 
methanol  are  usually  equal  to  or  less  than  unity.  Therefore, 
raa  and  mc  should  be  equal  or  less  than  ba  and  bc ,  respec¬ 
tively.  However,  the  mass-transfer  rates  of  methanol  in  the 
anode  and  oxygen  in  the  cathode  would  be  affected  by  other 
processes  such  as  water  flooding  on  the  cathode,  blocking 
of  the  methanol  pathway  by  CO2  bubbles  on  the  anode. 
Therefore,  ma  and  mc  in  Table  4  and  Table  5  obtained  by 
curve  fitting  of  the  experimental  data  were  different  from 
their  theoretical  value. 

(iv)  n a  and  nc  should  be  equal  to  the  reciprocals  of  the 
anodic  and  cathodic  limiting  current  density  and  their 


values  should  decrease  with  increasing  temperature.  Our 
results  in  both  Tables  4  and  5  are  consistent  with  this 
prediction. 

(v)  Parameters  related  to  methanol  crossover  0*rOssOCP’ 
Across, a,  and  Across, c)  did  not  appear  in  any  semi-empirical 
equation  proposed  in  the  literature  before.  /*ross  OCP  is  the 
equivalent  current  density  of  the  methanol  crossover  flux 
from  the  anode  at  open  circuit.  The  value  of  ^rossOCP 
should  decrease  as  the  temperature  is  lowered  or  the 
methanol  solution  is  diluted  [18,19].  As  for  £cr0Ss,a,  its 
value  is  dependent  on  the  methanol  crossover  flux.  With 
a  given  anodic  limiting  current,  larger  /*ross  OCp  will  raise 
Across, a  and  reduce  the  cell  voltage.  In  contrast,  £Cross, c  rep¬ 
resents  the  intrinsic  voltage  loss  at  cathode  due  to  methanol 
crossover.  If  0  is  a  characteristic  value  of  the  used  catalyst, 
Across, c  could  be  regarded  as  a  constant  for  a  given  ME  A. 
For  this  reason,  £cross? c  should  be  independent  of  /*ross  OCP 
and  our  fitting  results  (Tables  4  and  5)  do  demonstrate  that 
Across, c  is  nearly  constant  at  40  and  60  °C  for  the  same 
MEA. 

4.2.  Cathodic  and  anodic  polarization  curve 

Fig.  6  shows  the  polarizations  of  the  cathode  and  anode  indi¬ 
vidually  based  on  our  model.  The  anodic  polarization  curve  was 
sensitive  to  temperature  change,  which  implies  that  by  raising 
temperature  the  improved  power  density  of  the  whole  cell  is 
mainly  due  to  the  anodic  side’s  improvement.  By  analyzing  the 
anodic  parameters,  we  believe  the  anodic  polarization  is  domi¬ 
nated  by  the  changes  of  £o,a  and  Across, a  responding  to  different 
temperatures.  In  addition,  the  cathodic  polarization  curve  shown 
in  Fig.  6  includes  methanol  crossover,  in  contrast  to  the  cathodic 
polarization  shown  in  Fig.  3,  which  does  not  include  methanol 
crossover.  This  distinct  difference  between  these  two  cathodic 
polarization  curves  represents  the  crossover  effect  significantly. 
Apparently,  £cr0Ss, c  and  /*ross  OCP  are  the  key  parameters  of  the 
cathodic  performance  curves.  Besides,  /*ross  OCP  which  is  sensi¬ 
tive  to  temperature  is  also  tightly  associated  with  Across, a-  Again, 


Fig.  6.  The  individual  polarization  behavior  of  anode  and  cathode  at  different 
temperatures  (40  and  60  °C). 
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this  proves  that  the  anodic  and  cathodic  crossover  effects  are 
interrelated. 

4.3.  Cell  voltage  loss  due  to  individual  overpotentials 

Fig.  7  is  a  breakdown  of  overpotentials  in  our  ME  A  experi¬ 
ments.  It  clearly  illustrates  the  relative  contribution  of  individual 
overpotentials.  The  activation  overpotentials  of  both  electrodes 
are  the  major  cause  of  voltage  loss,  implying  the  importance 
of  catalyst  performance.  The  cathodic  concentration  overpoten¬ 
tial  is  smaller  than  the  concentration  overpotential  on  the  anode. 
This  is  due  to  excessive  air  flow  on  the  cathode.  The  mass- 
transfer  effect  has  a  minor  effect  on  the  overpotential  loss.  In 
Fig.  7,  the  voltage  loss  due  to  the  internal  resistance  becomes 
significant  in  the  high-current  region.  In  our  model,  the  total 
voltage  loss  caused  by  methanol  crossover  can  be  defined  as  the 
difference  between  potential  calculated  from  Eq.  (14)  and  that 
calculated  from  Eq.  (15)  which  eliminates  the  crossover  effect 
(<5  =  0,  /cross, ocp  =  0?  Across, a  =  0  Across, c  =  0)-  Besides,  the  anodic 
and  cathodic  crossover  effects  can  be  further  separated  by  com¬ 
paring  the  anodic  or  cathodic  part  between  Eqs.  (14)  and  (15). 

/^cell  =  Ec  Ea  ^celfl 

=  E0,c  -  bc  ln(/*)  +  mc  ln(l  ncC)  E0 ,a 

-  K  ln(/*)  +  ma  ln(l  -  naC)  -  Rce ui*  (15) 


4.4.  Impact  of  individual  parameter  on  the  cell 
performance 

Activation  overpotential,  internal  resistance,  concentration 
overpotential  of  oxygen  and  methanol,  and  methanol  crossover 
are  the  main  causes  for  cell  voltage  loss.  With  proper  ME  A 
design  and  fabrication  process,  these  overpotentials  could  be 
reduced.  Figs.  8-12  illustrate  how  individual  parameters  can 
affect  the  cell  performance  (output  voltage,  ECQ n,  and  power  den¬ 
sity,  P ).  For  simplicity,  Eq.  (15)  was  used  and  only  changes  in 
the  cathodic  parameters  were  introduced,  and  the  anodic  per¬ 
formance  was  calculated  by  fixing  the  anodic  parameters  at 
£o,a  =  0-3  V,Z?a  =  0.02  V  e_1,raa  =  0.1  V  e_1  and  na  =  2  cm2  A-1. 
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Fig.  8.  The  effects  of  £o,c  (^c  =  0.02  Ve  1,mc  =  0.05Ve  1,^c  =  lcm2A  1  and 
RCe\\  =  0.5  Q,  cm2). 


Fig.  9.  The  effects  of  bc  (%  =  0.8V,  mc  =  0.05Ve  1 ,  nc  =  lcm2A  1  and 
Rce\\  =  0.5  Q,  cm2). 


In  fact,  changes  of  anodic  parameters  would  yield  similar  effects 
on  cell  performance.  Both  cell  voltage  ( E )  and  power  density  (P) 
were  plotted  as  functions  of  current  density  (i  )  on  the  same  dia¬ 
gram. 


Fig.  7.  The  ratio  of  every  overpotential  in  the  entire  current  range  at  60  °C. 


Fig.  10.  The  effects  of  Rce\\  (£o,c  =  0.8  V,  bc  =  0.02  V  e  1 ,  mc  =  0.05  V  e  1  and 
nc  =  1  cm2  A-1). 
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Fig.  11.  The  effects  of  mc  (£o,c=0.8V,  bc=  0.02  Ve  1 ,  nc  =  lem2A  1  and 
Rce\\  =  0-5  £2  cm2). 


Fig.  13.  The  effects  of  £Cross, c  (%  =  0.8V,  £>c  =  0.02  V e  1 ,  rac=0.05Ve  1 , 
Rce\\  =  0-5  92  cm2,  nc  =  1  cm2  A-1  and  z'*ross  0cp  =  0.05  cm2  A-1  (the  calculated 
Across, a  =  0.044  V)). 


Besides  operating  condition,  parameter  £o,c  and  bc  were 
heavily  related  to  the  catalyst’s  inherent  properties  including 
its  composition,  structure  and  particle  size.  As  shown  in  Fig.  8, 
variation  of  £o,c  could  greatly  change  the  cell  voltage  at  a  given 
current  density.  The  effect  of  bc  on  the  cell  voltage  (Fig.  9) 
would  be  more  significant  in  the  low  current  density  region  and 
less  in  the  high  current  density  region.  By  comparing  Fig.  8  with 
Fig.  9,  £o,c  was  more  influential  on  the  power  density  than  bc. 
So  the  increase  of  £o,c  should  be  the  main  criterion  for  catalyst 
improvement. 

Fig.  10  illustrates  the  effect  of  internal  resistance  on  cell  volt¬ 
age  and  power  density.  Since  the  voltage  loss  due  to  the  internal 
resistance  is  proportional  to  the  current  discharge,  the  voltage 
loss  is  significant  in  the  high  current  density  region.  In  Fig.  10, 
the  current  density  where  the  cell  delivered  maximum  power 
density  was  shifted  to  higher  current  densities  as  the  internal 
resistant  decreased. 

The  porosity,  hydrophobicity  and  the  thickness  of  the  gas 
diffusion  layer  and  catalyst  layer  all  affect  the  transfer  rate  of 
the  reacting  species  as  well  as  the  value  of  mc  and  nc  as  shown 
in  Fig.  11  and  Fig.  12.  Water  flooding  on  the  cathode  causes 
increases  in  mc  and  nc .  Reducing  mc  and  nc  tend  to  increase  the 


Fig.  12.  The  effects  of  nc  (£o,c  =  0.8V,  &c  =  0.02Ve  1 ,  mc  =  0.05Ve  1  and 
Rce\\  =  0-5  92  cm2). 


maximum  power  density  and  makes  it  shift  to  the  high-current- 
density  region.  However,  it  seems  the  effects  of  mc  and  nc  are 
minor  because  the  anodic  limiting  current  is  small  and  dominates 
the  overall  mass-transfer  overpotential. 

*  cross, OCP’  'E’cross.a,  and  £Cross,c  are  related  to  methanol 
crossover.  Their  influences  on  cell  voltage  and  power  density 
are  shown  in  Figs.  13  and  14,  respectively.  In  Fig.  13,  £cr0Ss, c  is 
similar  to  the  additional  activation  overpotential  at  the  cathode. 
As  its  definition,  £Cross, c  is  strongly  associated  with  bc  and  0. 
Accordingly,  a  more  methanol-tolerant  catalyst  may  be  helpful 
to  improve  the  cell  performance  because  a  larger  0  could  reduce 
Across, c-  Besides,  the  increase  of  bc  is  also  beneficial  for  the 
diminution  of  £Cross, c-  On  the  other  hand,  /*ross  OCp  and  Across, a 
are  reciprocal  parameters.  Fig.  14  shows  the  cell  performances 
with  different  values  of  /*ross  OCP  and  their  corresponding  Across, a 
with  unchanged  limiting  current  density  of  the  anode.  Moreover, 
smaller  raa,  which  means  better  mass  transfer  at  the  anode,  will 
be  helpful  to  mitigate  crossover  effect. 


Ecell  P 


Fig.  14.  The  effects  of  z*ross0 CP  and  £Cross,a  (%  =  0.8V,  bc=  0.02 Ve  \ 
mc  =  0.05  Ve-1,  RCq\\  =  0.5  92  cm2,  nc  =  1  cm2  A-1  and  Across, c  =  —0.05  V). 
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5.  Conclusion 

A  useful  semi-empirical  model  was  established  for  a  DMFC 
in  this  work.  This  model  can  analyze  anodic  and  cathodic  voltage 
losses  individually  so  that  the  controlling  factors  can  be  iden¬ 
tified  for  a  specific  system.  Parameters  involved  in  this  model 
include: 

(i)  Overpotentials  of  methanol  oxidation  reaction  on  the  anode 
(Eo, a,  and  hA). 

(ii)  Overpotential  of  oxygen  reduction  reaction  on  the  cathode 
(£b,o  and  bc). 

(iii)  Overpotential  losses  due  to  internal  resistance  (RCell)- 

(iv)  The  effects  of  methanol  crossover  on  the  both  anodic  and 
cathodic  potentials  (/*ross  OCP’  Across, a,  and  £Cross, c). 

(v)  Overpotentials  due  to  the  depletion  of  reacting  species,  such 
as  oxygen  ( nc  and  mc )  and  methanol  (na  and  ma). 

The  impact  and  physical  significance  of  major  parameters 
can  be  summarized  as  follows: 

(i)  In  the  low  current  density  region,  the  variations  of  EqiC  and 
Across, c  are  the  major  voltage  losses  at  the  cathode,  and  £o,a 
and  Across, a  play  the  same  role  as  £o,c  and  £Cross, c  at  the 
anode. 

(ii)  In  the  intermediate  current  density  region,  ohmic  and  con¬ 
centration  overpotential  are  the  key  reasons  causing  the 
voltage  drop. 

(iii)  With  increases  in  the  value  of  R ,  m  and  n ,  the  maximum 
power  density  decreases  and  the  current  density  accordingly 
drops. 

(iv)  /*ross  OCP  and  Across, a  are  inherently  coupled.  In  other  words, 
the  anodic  voltage  is  deeply  influenced  by  the  amount  of 
methanol  crossover  flux. 

(v)  Improvement  of  both  catalysts  and  mass  transfer  would  be 
useful  to  reduce  impact  of  the  crossover  effect,  as  charac- 
terized  by  Coss.OCP-  Across, c  and  £cross, a  in  Eq.  (14). 

The  performance  of  a  single  cell  can  be  analyzed  with  this 
model  to  determine  the  factors  that  cause  major  voltage  losses 
under  various  operating  conditions,  so  that  the  performance  of 
a  DMFC  can  be  optimized  without  too  much  trial-and-error. 
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